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ABSTRACT 

Context. The mass domain where massive extrasolar planets and brown dwarfs lay is still poorly understood. Indeed, not even a 
clear dividing line between massive planets and brown dwarfs has been established yet. This is partly due to the paucity of this kind 
of objects orbiting close to solar-type stars, the so-called brown dwarf desert, that hinders setting up a strong observational base to 
compare to models and theories of formation and evolution. 

Aims. We search to increase the current sample of massive sub-stellar objects with precise orbital parameters, and to constrain the 
true mass of detected sub-stellar candidates. 

Methods. The initial identification of sub-stellar candidates is done using precise radial velocity measurements obtained with the 
SOPHIE spectrograph at the 1.93-m telescope of the Haute-Provence Observatory. Subsequent characterisation of these candidates, 
with the principal aim of identifying stellar companions in low-inclination orbits, is done by means of different spectroscopic diag- 
nostics, as the measurement of the bisector velocity span and the study of the correlation mask effect. With this objective, we also 
employed astrometric data from the Hipparcos mission and a novel method of simulating stellar cross-correlation functions. 
Results. Seven new objects with minimum masses between ~ 10 Mjup and ~ 90 Mj^p are detected. Out of these, two are identified as 
low-mass stars in low-inclination orbits, and two others have masses below the theoretical deuterium-burning limit, and are therefore 
planetary candidates. The remaining three are brown dwarf candidates; the current upper limits for their the masses do not allow 
us to conclude on their nature. Additionally, we have improved on the parameters of an already-known brown dwarf (HD137510b), 
confirmed by astrometry. 

Key words, planetary systems ~ techniques: radial velocities - stars: brown dwarfs - stars: individual: HD14651, HD16702, 
HD22781, HD30246, HD9230, HD137510, HD156279, HD167215 



1. Introduction 

More than 15 years after the discovery of the first extrasolar 
planets and brown dwarfs (see, for example, [Mayor & Queloz 
[1995 Rebolo et al.|1995 i, a number of issues remains open re- 
garding the relation between these two types of objects. Principal 
among them is the problem of distinguishing members of one 
class from members of the other. As discussed by, for exam- 
ple, |Basri_&^rown (2006| the criteria used to classify these ob- 
jects can concern their physical characteristics, and their origin, 
among other things. 

The common idea that brown dwarfs are "failed" stars that 
only fusion deuterium, and that are not able to sustain a con- 
stant luminosity throughout their lifetimes (e.g. Burrows et al. 
1997 Baraffe et al. 2003 i, while planets are objects that are sim- 



* Based on observations collected with the SOPHIE spectrograph 
on the 1.93-m telescope at Observatoire de Haute-Provence (CNRS), 
France, by the SOPHIE Consortium (program 07A.PNPCONS). 



ply too light to fusion deuterium seems to imply that a defini- 
tion based on physical characteristics would be the best adapted 
one. Traditionally, the deuterium-burning mass was taken to be 
13 Mjup ( [Burrows et al.|200i ), but there exists an important un- 
certainty related to the equation of state used in these computa- 



tions. Additionally, more recent models by Spiegel et al. (201 1 
show that the deuterium-burning mass limit depends on the he- 
lium abundance, the initial deuterium abundance, and the metal- 
licity of the model, among other things. Therefore, the authors 
find that the deuterium-burning mass can range from 1 1 Mjup 
up to 16 Mjup. Choosing a single mass as the limit between 
both categories may therefore lead to an oversimplification of the 
scenario. Additionally, since the minimum mass needed to start 
gravitational collapse, either in a disk or in a molecular cloud, is 
probably around a few Jupiter masses (Whitworth & Stamatellos] 
2006), and since the core-accretion models predict objects with 
masses as high as a few tens of Jupiter masses ( jMordasini et al.| 
|2009 1, there surely exists an important overlap in mass for ob- 
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jects having different origins. For example, |Kirkpatrick et al. | 
( [2006 ') have reported the discovery of a L-type brown dwarf with 
a mass that could be as small as 6 Mjup. This might make the dis- 
tinction by mass not completely satisfactory. 

On the other hand, a criterion based on formation history 
would stand on the fact that brown dwarfs are thought to form 
as stars (see the reviews by Luhman et al. 2007 'Whitworth et al. 



2007 ), i.e. by gravitational collapse and fragmentation of molec- 



ular clouds on dynamical timescales, while planets are widely 
believed to form from accretion of a solid core in a protoplan- 
etary disk, a process believed to be much slower ( [Pollack et al 



|1996 ; Alibert et al. 2005 ). This could in principle allow us to dis- 
tinguish ones from the others, but this is not without a few draw- 
backs. First of all, it is practically impossible to trace unambigu- 
ously the history of an object back to its origins. Additionally, 
there is not a complete consensus about the mechanisms rul- 
ing formation processes, neither for brown dwarfs, nor for mas- 
sive planets. For example, massive planets are also proposed to 
be able to form by disk collapse ( |Boss|1997 1 or by protostellar 
cloud fragmentation followed inward migration by disk capture 
( |Font-Ribera et al. 2009). A definition based on formation his- 
tory might then prove provisional. 

To produce a consistent picture of brown dwarfs and extra- 
solar planets, a solid observational base for formation and evo- 
lution models is crucial. Constraining the multiplicity properties 
of these objects, such as frequency, separation, and mass ratio 
distributions, should permit distinguishing between different for- 
mation models. 

In this regard, most of the known extrasolar planets have 
been discovered by high-precision radial velocity (RV) surveys 
(see exoplanet.eu, and Schneider et al.,,2011, ). Due to the large 
reflex motion they induce in the stars they monitor, these sur- 
veys should easily discover brown dwarfs companions and prop- 
erly characterise their orbital elements, and frequency. However, 



the brown dwarf desert ( [Marcy & Butler] |2000| |Grether & 



Lineweaver||2006 1, i.e. the paucity of massive sub-stellar com- 
panions to solar-type stars in close orbits, has hindered this en- 
deavour, resulting in only few tens of such objects known to date 
( [Sozzetti &"D "esidera 20101 |Sahlmann et al. 2011b). Indeed, the 
frequency of close-in brown dwarf companions around sun-like 



stars has been recently determined to be below 0.6% ( Sahlmann 
let al. 2011b| l. Additionally, the RV method provides a lower limit 



for the mass only, which means that a fraction of the brown 
dwarf candidates detected by RV surveys are actually low-mass 
stars in low-inclination orbits. As a consequence, studies of mas- 
sive planets and brown dwarf properties are faced with a small 
and uncertain sample of objects, and therefore all efforts to ei- 
ther increase the number of these objects or to determine their 
true mass accurately should prove valuable. 

Different authors have attempted to measure the true mass 



of RV candidates using astrometric measurements (Halbwachs 
eTaLllSOOO, Zucker & Mazeh 2001, Sahlma mTet al.||201 Ibj 



Reffert & Quirrenbach 2011). Although most of these analyses 



fail to yield direct measurements of the mass, and are limited by 
the precision of their measurements to upper limits to the mass, 
they have proven efficient in discarding relatively massive stellar 
companions in low-inclination orbits. 

In this paper, we report the detection of 7 new objects 
with minimum masses between 10 and 90 Mjup, orbiting main- 
sequence stars of spectral types F, G and K. The detections is 
based on precise radial velocity measurements obtained as part 
of the sub-progamme 2 of the SOPHIE search for northern extra- 
solar planets (Bo uchy et al.|2009a ). We have already announced 
three of these objects, albeit with preliminary parameters (Diaz 



|et al.|20l"T| ). Additionally, we improve on the orbital parameters 
of an already known brown dwarf, HD 1375 10 b. Finally, our 
attempts to measure the true mass of these candidates are also 
described. 

In Sect. 12] we describe the observations and the data reduc- 
tion, in Sect.[3]we describe the methods used to analyse the data, 
fit the orbits, and estimate the stellar parameters of the primaries. 
We also describe the techniques employed to constrain the true 
nature of the companions, including the analysis of astrometric 
data from the Hipparcos mission. In Sect.|4]we report the results 
obtained for each object, which are roughly classified using the 
mass-based criterion described above into planetary candidates, 
brown dwarf candidates and stellar-mass objects. The results are 
then discussed in Sect. [5] Finally, in Sect. [6] we summarise the 
results and present our conclusions. 



2. SOPHIE Radial Velocities 

2.1. Observations 

Radial velocity measurements were performed with SOPHIE, 
a high-resolution echelle spectrograph fiber-fed from the 
Cassegrain focus of the I.93-m telescope at the Haute-Provence 
Observatory (OHP). It is installed in a temperature-stabilised 
environment and the dispersive elements are kept at constant 
pressure in order to provide high-precision radial velocities 
( jPerruchot et al.|2008 ). SOPHIE spectra cover most of the visi- 
ble wavelength range in 39 spectral orders. 

Observations were performed in high-resolution mode which 
produces spectra with resolving power of /l/A/l ~ 75 000 at 
550 nm. Two optical fibers are used in this mode: one of them 
records the stellar spectrum while the other can be used to si- 
multaneously calibrate the spectrum with a Thorium Argon lamp 
{thosimult mode) or to record the background sky spectrum dur- 
ing the observations (objAB mode). Most of the observations re- 
ported here were obtained in the objAB mode, since the stability 
of the instrument is high enough to allow for calibrations to be 
obtained every 2 hours without degrading the sought velocity 
precisiorQ The observations performed in thosimult mode date 
from the first season of observations with SOPHIE when its sta- 
bility had not yet been confirmed. Moreover in the objAB mode, 
the sky fiber can be used to correct for scattered light contami- 
nating the star fiber (e.g. [Barge et al.|2008[|Pollacco et al.|2008 
Santerne et al..2011) 



The targets reported in this paper are listed in Table [T] We 
include the date of the first and last observation of each target, 
the time span of observations and the number of measurements 
obtained. In the last column of Table[T]we indicate the Section in 
which the results of each target are presented. The exposure time 
was varied in order to reach a S/N of 50 per resolving element 
under varying weather conditions. The exposure time and S/N 
of each individual observation are reported in Table [2] available 
on-line. 



2.2. Data reduction 

The spectra were reduced and extracted using the SOPHIE 
pipeline (Bouchy et al. 2009a| ), which performs all the steps 
from the raw images to the wavelength-calibrated 2D spectra. 
These are then correlated using numerical binary masks corre- 



' The same is not true for other SOPHIE progammes, which require 
measurements with the highest possible precision (see [Bouchy et al^ 
2009a| . 
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Table 1. Target characteristics and summary of observation log. 



Target 


Coordinates 


m. 


Sp. Type 


Dist." 




Observations 




Sect. 


Name 


R.A. 

(hli:mm:ss) 


DEC 

(-:':") 






(pc) 


First Obs.* 


Last Obs.* 


Time Span 
(days) 


N 




HD 14651 


02:22:00.9 


+04:44:48.3 


8.3 


G6 


39.6 


2006-11-04 


2010-09-14 


1410.1 


15 


4.2 




HD 16702 


02:42:58.7 


+58:11:52.1 


8.3 


GO 


45.2 


2007-09-22 


2010-11-02 


1137.9 


22 


4.3 




HD22781 


03:40:49.5 


+31:49:34.7 


8.8 


Kl 


31.8 


2006-11-12 


2010-09-05 


1392.1 


33 


4.1 




HD30246 


04:46:30.4 


+ 15:28:19.4 


8.3 


Gl 


51.5 


2006-11-12 


2010-10-18 


1436.9 


25 


4.2 




HD92320 


10:40:56.9 


+59:20:33.0 


8.4 


G3 


41.5 


2006-12-21 


2010-11-27 


1437.1 


16 


4.2 




HD137510 


15:25:53.3 


+ 19:28:50.6 


6.3 


F9 


41.8 


2008-03-24 


2010-08-27 


886.7 


13 


4.2 




HD156279 


17:12:23.2 


+63:21:07.5 


8.1 


G6 


35.4 


2010-05-14 


2011-01-23 


254.2 


15 


4.1 




HD167215 


18:12:59.4 


+28:15:27.4 


8.1 


F7 


78.6 


2009-07-29 


2010-10-09 


437.0 


14 


4.3 





Notes. Distance from Hipparcos catalogue ( Ferryman et al.|l997| . Universal Time date at beginning of observation. 



sponding to different spectral types (G2, K5 and M4). The cross- 
correlation function (CCF) averaged over all 39 spectral orders is 
then fit with a gaussian function to obtain the RV, the full-width 
at half-maximum (FWHM), and contrast relative to the baseline. 
Additionally, the bisector velocity span of the CCF is computed 
in the way described bylQueloz et al. ( |2001) for further analysis, 
as well as the radial velocity computed considering only the blue 
or red halves of each spectral order, which we call RVuue, and 
RViei, respectively (see ^3.4.1 1. For reference, in Table[3]we re- 
port for all targets studied in this paper the numerical mask clos- 
est to the measured spectral type (^TTJ, and the median values 
of the CCFs parameters: RV, FWHM and Contrast. The errorbars 
reported correspond to the standard deviation of each magnitude. 

The complete radial velocity datasets obtained with the mask 
closest to the spectral type of each target are reported in Table |2] 
where the uncertainties include the photon noise error, estimated 
with the method by Bouchy et al.| (2001 1, the error in the wave- 
length calibration (~ 2 ms"'), and an additional systematic error 
of 3 m s which was added quadratically to the other sources of 
uncertainty. The systematic errors arise from imperfect guiding, 
the so-called seeing effect ( jBoisse et al.|[2011a , see also 5 3.4.1 
of the present paper), uncorrected velocity drifts of the spectro- 
graph, and other still unidentified sources. 

In April 2007, SOPHIE was upgraded with a bandpass filter 
to reduce the contamination of the stellar spectra by Argon lines 
coming from the ThAr lamp during observations in thosimult 
mode. For three stars of our sample some observations were per- 
formed in thosimult before the upgrade: HD14651, HD22781, 
and HD30246. These observations are identified in Table |2] and 
they exhibit a clear contamination in the red part of the spectra. 
To try to mitigate the effects of this contamination, the 6 reddest 
spectral orders were not considered for the computation of the 
average CCF and the radial velocity. To maintain the uniformity 
of the reduction, the remaining spectra of each star were also 
reduced in the same way. The resulting velocities are slightly 
different from those obtained using all spectral orders for these 
same spectra, but no systematic shift exists between both ob- 
servation modes. In general, the orbital fits performed using the 
corrected velocities exhibit a reduced scatter. The remaining ob- 
servations were performed in objAB mode. 



Table 3. Cross-correlation function parameters. Median values 
are reported; the error intervals correspond to the standard devi- 
ation of each magnitude. 



Target 


Mask 


RV 


FWHM 


Contrast 






[km s-'] 


[kms-'] 


[%] 


HD14651'' 


G2 


52.2 


7.53 ± 0.04 


43.1 ±0.6 


HD 16702 


G2 


5.5 


7.58 ± 0.02 


35.8 ±0.4 


HD22781 " 


K5 


8.2 


6.60 ± 0.03 


34.3 ± 1.7 


HD30246'' 


G2 


41.2 


9.34 ±0.14 


35.1 ±2.4 


HD92320 


G2 


1.1 


7.70 ± 0.02 


38.3 ±0.5 


HD137510 


G2 


-6.7 


12.14 ±0.06 


29.1 ±0.4 


HD156279 


G2 


-20.5 


7.64 ±0.01 


47.1 ±0.5 


HD167215 


G2 


-44.6 


9.74 ± 0.04 


22.2 ±0.7 



Notes. Only spectra acquired in objAB mode are considered for the 
contrast and the FWHM. 



tational velocity (y sin/, where / is the inclination of the ro- 
tation axis of the star with respect to the line of sight) using 
the calibrations described in detail by |Boisse et al. ( 2010| l. The 

on 



estimated uncertainty due to these calibrations is 1 kms"^ 
V sin / and about 0.09 dex on [Fe/H]. Additionally, a more de- 
tailed spectral analysis was carried out using the SOPHIE spec- 
tra obtained in objAB mode, and not contaminated by moonlight. 
The method ( Santos et al.||2004) l provides the effective temper- 
ature Teff, the surface gravity log g and the metallicity [Fe/H] 
of the observed stars. To obtain the stellar masses, we interpo- 
lated the isochrones by Girardi et al. (2000i, using the method 
described by |da Silva et al. |2006) l. The values derived from this 
analysis are reported in Table |4| The uncertainties for the stellar 
masses range from 2.6% to 5.8%. However, these are the formal 
uncertainties, which do not take into account the systematic un- 
certainties in the models, and we therefore decided to adopt a 
conservative uncertainty of 10% on the stellar mass. The spec- 
tral types reported in Table [T] are obtained from interpolating the 
tables from |Cox| (|2000 ) to the values of T^.^ obtained. We note 
that in all cases, the spectroscopic determination of the metallic- 
ity agrees with the method based on the CCF ( [Boisse et al. 2010| l 
to within 1-cr. 



3. Data Analysis 

3.1. Stellar Parameters 

The CCF obtained ( |2.1| i has been used to obtain an estimate 
of the stellar metallicity ([Fe/H]) and the projected stellar ro- 



3.2. Stellar activiy 



We computed the log/?™ activity index (Noyes et al. 



1984 1 



on each spectrum of each target using the rnethod described by 
[Boisse et al.| ( |2010[ ). In Table |4] we report the median value of 
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Table 4. Stellar parameters of target stars. 



Parameters 


HD 14651 


HD 16702 


HD22781 


HD30246 


Effective temperature [K] 


5491 ± 26 


5908 ± 25 


5027 ± 50 


5833 ± 44 


vs'ml [kms"']" 


3.07 


2.56 


1.73 


4.96 


[Fe/H] 


-0.04 ± 0.06 


-0.12 ± 0.06 


-0.37 ± 0.12 


0.17 ± 0.10 


log g 


4.45 ± 0.03 


4.46 ± 0.03 


4.60 ± 0.02 


4.39 ± 0.04 


Ms [Mo]* 


0.96 ± 0.03 


0.98 ± 0.04 


0.75 ± 0.03 


1.05 ±0.04 


Age [Gyr] 


[0.2 - 6.6] 


[0.1-4.6] 


[0.6 - 8.9] 


[0.2 - 3.7] 


Age;;;^^ [Gyr]'^ 


> 5.0 


[0.6- 1.1] 


> 4.6 


[0.4 - 0.9] 




-5.02 ±0.11 


-4.54 ± 0.05 


-5.02 ±0.14 


-4.50 ± 0.06 


Parameters 


HD92320 


HD137510 


HD 156279 


HD167215 


Effective temperature [K] 


5664 ± 24 


6131 ±50 


5453 ± 40 


6351 ± 68 


vs'ml [kms"']" 


3.13 


7.42 


2.51 


4.94 


[Fe/H] 


-0.10 ±0.06 


0.38 ±0.13 


0.14 ±0.01 


-0.21 ±0.14 


log g 


4.48 ± 0.03 


4.02 ± 0.04 


4.46 ± 0.03 


4.19 ±0.07 


Ms [M^f 


0.92 ± 0.04 


1.36 ±0.04 


0.93 ± 0.04 


1.15 ±0.07 


Age [Gyr] 


[0.2 - 6.6] 


[2.3 - 3.0] 


[0.4 - 7.7] 


[2.0 - 5.2] 


Age«^_^ [Gyr]^ 


[2.2 - 4.7] 


> 5.4 


> 5.4 


> 8.0 


log^HK 


-4.80 ±0.08 


-5.09 ±0.18 


-5.01 ±0.06 


-5.25 ±0.19 



Notes. Estimated from CCF as described in Sect. 3. 1 Uncertainty is estimated to about 1 km s"' in v sin / (see also Boisse et al. 2010 1. Formal 
errors issued from isochrone fitting are reported. A more realistic 10% was used for the rest of the analysis. ' Age estimated from the age-activity 
relation by |Mamajek & Hillenbrand| ( |2008^ . 



the index for each target and its standard deviatiorj^ computed 
using the spectra obtained in objAB mode only, since the con- 
tamination by the ThAr lamp present in the rest of the spectra is 
expected to change the depth of the spectral lines. The obtained 
values were used to estimate the age of the host stars using the 
relation by M amajek & Hillenbrand| ( 2008 1. We report the results 
in Table [4] The activity-determined ages are in agreement with 
the values obtained from isochrone fitting, except for HD 1 375 1 
and HD 1672 15. These stars are the two least active stars in the 
sample, and their log/?^j^ are either at the limit of or outside 



the range of applicability of the IVlamajek & Hillenbrand ( 2008| l 



formula, which is valid only for log/?^^ > -5.1. 

It is known that stellar activity is a source of uncertainty 
in the RV measurements (e.g. ]Queloz et al.||2001l |Boisse et al. 
|2011b| ). In some cases (Boiss^et al. 2010; Melo et al. 20()7|l, the 
RV jitter produced by stellar activity has been corrected, with 
the result of improved determinations of the orbital and plane- 
tary parameters. This is done by decorrelating the residuals of 
the orbital fit with the bisector velocity span. For all our sample 
stars, we checked if the residuals of the best-fit Keplerian or- 
bit exhibited any dependence with the measured bisector span. 
In no case a correlation is found with a significance above 2 
cr. In the case of HD30246, the most active star in our sample 
(log/?Jjj^ - -4.5), a slight hint of negative correlation is seen, 
with a fit slope ofb = -0.13 + 0.08. On the other hand, HD16702 
has a similar activity level (log R'^^^ - -4.54), but the residuals 
do not exhibit any correlation with activity index RUur- However, 



as will be discussed in ^4.3[ the companion to HD 16702 is most 
likely a low-mass star, which might contaminate the determina- 
tion of R[if^- We note that the stars studied by Boisse et al. ( 2010 1 
(HD189733) andl JVIelo et aL] ( p007) (HD102ly5) exhibit activity 
levels similar to that of HD30246 and HD 1 6702, but their data is 
much more precise, which might permit detecting subtler trends. 



3.3. Moon contamination 

Whenever present, moonlight scattered in the Earth's atmo- 
sphere can pollute the observed stellar spectrum. If the Moon's 
RV is close to the velocity of the observed star, moonlight can 
bias the measurement of the velocity. Prior to performing a kep- 
lerian fit to the data, we have verified that none of the RV mea- 
surements were biased by this effect. For an observation to be 
polluted, two conditions must be fulfilled: in the first place, the 
scattered moonlight must contribute significantly to the light en- 
tering the spectrograph's fiber. For the observations in objAB 
mode, this can be addressed using fiber B. We considered con- 
servatively that the moonlight level is enough to pollute a stel- 
lar spectrum if both the S/N in fiber B was larger than 3 and 
the contrast of the corresponding CCF was larger than 1%, i.e. 
if the CCF peak of the Moon is detected in the sky fiber, even 
if marginally. The second condition is that the Moon RV must 
be close to the stellar velocity. To study this, we considered the 
Baricentric Earth Radial Velocity (BERV) on the day of the ob- 
servation, and in the direction of the target, which is within ~ 1 
km s ' of the Moon Radial Velocity. The Moon RV is considered 
to be close to the measured RV if the BERV was less than ~ Icr 
away from the measured stellar peak, where cr is the width of 
the the CCF. For those cases in which the second fiber was dedi- 
cated to obtain the spectrum of a ThAr lamp simultaneously, for 
which no sky monitoring is possible, the presence of the moon 
was assumed whenever its illuminated fraction was larger than 
40%. Besides this, the same criterion based on the BERV was 
applied. 

With these criteria, none of the reported velocities is pol- 
luted, due mainly to the fact that the targets are bright. However, 
for limiting cases, to confirm that the effect of moonlight pol- 
lution was negligible, we performed the Moon correction as is 
described by [Barge et aE] ([2008]) andlPollacco et al.| ( |2008) l for 
SOPHIE and by |Bonomo et al.[ ( |201o| for HARPSTin all cases. 



^ The standard deviation of log R'^^, was computed using the median 
absolute deviation, which is less sensitive to the presence of outliers. 



we found that the corrected velocities were in excellent agree- 
ment with the uncorrected measurements, which were therefore 
retained. 
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3.4. Keplerian Orbits 



3.5. Resolving fhe nature of tine companions 



The radial velocity measurements of all the targets exhibit vari- 
ations greatly in excess of the measurement precision. The RVs 
were fit using a single Keplerian orbit model. A genetic al- 
gorithm was employed to pick up the starting point for the 
Levenberg-Marquardt method. The best-fit parameters are re- 
ported in Table [6j together with the 68% confidence interval, 
obtained from 5000 Monte Carlo simulations. In most cases, the 
histograms of the obtained parameters have a single peak, and 
gaussian-like appearance, and therefore the errors reported can 
be taken as the l-cr intervals. In a few cases, however, double- 
peaked or highly asymmetric distributions were found. These 
cases are studied in detail for each object in ^ In Figure [T] we 
plot the RV curves for all stars observed together with the best- 
fit model and the residuals, after correction of the seeing effect 
(see below). 



3.4.1 . Correction of the seeing effect 

Radial velocities obtained with SOPHIE exhibit a systematic ef- 
fect thought to be caused by an incomplete radial scrambling in 
the optical fiber coupled with aberrations in the optical system, 
which lead to a variation in the illumination of the spectrograph's 
pupil. The consequence of this effect is an artificial variation in 
the measured velocity when the point spread function changes, 
as a consequence of seeing variations, focus variations, etc. This 
effect induces RV variations with an amplitude of about 20 m s"' 
and constituted the main limitation on the spectrograph's preci- 
sion until very recentljj^ 

Boisse et al. (201 la!) present a detailed description of this ef- 



fect and provide a way to diagnose and correct it. We follow the 
method described by these authors, but using a different variable 
to decoiTelate the data. Since simulations of the optical path of 
SOPHIE have shown that variations in the pupil's illumination 
produce a differential effect between the red and blue parts of 
each order ( Boisse etal.|2011al ), we have used the difference be- 
tween the velocities measured in each half of the CCD (see |2.1| i 
as an indicator of this effect. We define the pupil's illumination 
proxy 6rv as 

6rv - RVhiue - RVi-ed , 

where ^Vbiue and RVi-ed indicate the velocities measured using 
the blue and red part of the orders, respectively. 

Not all stars in our sample exhibit a clear correlation between 
the residuals and dny- This is also the case when the estimator 
from [Boisse et al.| ( |2011a| ) is used, but the reasons for this dif- 
ferential behaviour are currently under study. For the moment, 
we have decided to perform the correction only on the stars for 
which the Pearson's correlation coefficient between 6rv and the 
fit residuals is greater than 0.5. In this case the residuals are 
decorrelated and added back to the fit curve. The resulting data 
set is then used to find a new solution to the 1 -Keplerian fit, 
which is shown in Figure [T] The data of only three stars fulfil 
this criterion: HD16702, HD92320 and HD137510, for which 
the correlations between residuals and 6rv are shown in Fig. |6] 
For the rest of the sample, the coiTelation is not significant, and 
the data reported do not include any correction. 



The high minimum mass of most of the reported companions 
imply a non-negligible probability that the objects are actually 
low-mass stars in low-inclination orbits. To narrow down this 
possibility and, when possible, conclude on the nature of the de- 
tected objects, we have carried out a series of tests and simula- 
tions using the available RV data and astrometric observations 
from the Hipparcos satellite ( Ferryman et al.|1997| l. 



^ Tfiis has been solved by a recent upgrade of the fiber link ( Perruchot 
et al.|2011ll. 



3.5.1 . Bisector Analysis and CCF Simulations 

If the companion object is relatively luminous, we expect it to 
pollute the peak of primary star in the CCF, and to produce a 
variable distortion of its shape, which results in a variation of 
its bisector velocity span throughout the object's orbit. We have 
searched for significant variations in the bisector velocity span, 
as well as correlations between the bisector velocity span and 
RV measurements. The detection of such correlation indicates 
that the velocity signal is the product -at least in part- of a de- 
formation of the stellar CCF. This can either mean that the ve- 
locity signal is in fact a product of stellar activity -which would 
be unprecedented for the large RV amplitudes observed in these 
targets, and that could easily be identified using other diagnos- 
tics such as the correlation between RV and the log R'^^, index-, 
or that the companion is contributing non-negligibly to the ob- 
served CCF. For all the stars of our sample, except for HD 16702 
(see ^below), the bisector velocity span observations have neg- 
ligible variability and do not exhibit any significant correlation 
with radial velocity measurements. 

In the context of this paper, the lack of correlation between 
the bisector velocity span and the radial velocities can be used 
to constrain the mass of the companion by means of CCF sim- 
ulations. To exploit this possibility, we assume that both the 
host and the companion objects are dwarf stars, and that each 
contributes its own peak to the observed CCF. We model each 
CCF as a gaussian function. The parameters of the primary CCF 
set to the median values of the observed quantities. The area 
and FWHM of the secondary peak are varied according to the 
modeled mass ratio and rotational velocity, assuming the same 
metallicity as the primary star, and using the calibrations from 
Boisse et al. (2010 1. For each observed orbital phase, the pri- 
mary and secondary CCF are simulated at the expected position 
in RV space, and are added after normalisation using the flux 
ratio of the stars. Here we assume that the measured RV coiTe- 
sponds approximately well to the RV of the primary star, which 
holds for small mass ratios. For each simulated orbital phase, 
the RV, bisector velocity span, FWHM and contrast of the joint 
CCF are computed similarly to how is done for the real data, 
and the results are compared to the observations. The value of 
the;\f^ statistics is computed over a grid of masses and rotational 
broadening values for the secondary, and confidence regions are 
obtained from intersecting the surface with constant planes. 
In general, we have found that the bisector velocity span pro- 
vides strong constrains on the mass of the secondary, resulting 
in many cases in lower upper limits than those issued from as- 
trometric measurements. The results from these simulations are 
reported for each target in ^ 

3.5.2. Mask Effect 

Similar to the bisector analysis, the dependence of the amplitude 
of the RV variations with the numerical mask used to measure 
them can reveal the contribution of a secondary stellar object to 
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Table 5. Parameters of the Hipparcos astrometric observations 



Target 


HIP 


Sol. 


A'o.-b" 




A'Hip" 






type" 




(mas) 




HD 14651 


11028 


5 


13.4 


3.4 


59 


HD 16702 


12685 


5 


15.0 


5.3 


140 


HD22781 


17187 


5 


1.7 


4.0 


69 


HD30246 


22203 


5 


0.8 


3.5 


56 


HD92320 


52278 


5 


8.1 


6.4 


233 


HD137510 


75535 


5 


1.4 


1.9 


101 


HD 156279 


84171 


5 


8.7 


3.5 


110 


HD167215 


89270 


7 


1.8 


5.3 


144 



Notes. '"' '5': standard five-parameter astrometric solution; '7': solu- 
tion including first derivatives of the proper motion. Portion of the 
orbit covered by the Hipparcos measurements. Median precision. 
Number of observations. 
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Fig. 2. Residuals of a single-Keplerian fit to the radial velocities 
of HD22781, not including any additional drift. 



the flux entering the spectrograph ( Bouchy et al.|2009b Santos 
let al.|2002| i. 

When analysing our data in search of this effect, we proceed 
in the following way: first, we compute the best Keplerian fit 
as described in 5 3.4 using the best-matching numerical mask 
(Table |3]l. Then we adjust the RV measurements obtained with 
other masks, fixing all orbital parameters, except for the semi- 
amplitude K and the center-of-mass velocity y, which is ex- 
pected to change from mask to mask. Since the fit is linear on 
these two parameters, the errors are obtained directly from the 
covariance matrix. For all our targets, except for HD 16702 (see 



5 4.3 1, the amplitudes obtained in all three masks are compatible. 



sion of cTj^. Outliers in the IAD had to be removed because they 
can substantially alter the outcome of the astrometric analysis. 
We eliminated one outlier for HD30246 and HD92320 and two 
outliers for HD156279. 

An orbit signature is detected for HD16702, with a 2.5-cr 
significance, which provides a measurement of the real mass of 
the companions and places it in the stellar-mass range (see ^4.3| l. 
For the rest of the targets, no significant orbit is detected, but 
upper mass limits are derived. These are reported on a target 
basis in ^ The only exception is HD30246, for which not even 
an upper limit can be reliably derived from Hipparcos data (see 



3.5.3. Analysis of the Hipparcos astrometry 

All stars listed in T able [T] were cata logued by the Hipparcos 
astrometry satellite ( |Perryman et al. \991). We used the new 
Hipparcos reduction ( jF. van Leeuwen 2007 1 to search for sig- 
natures of orbital motion in the Intermediate Astrometric Data 



(IAD). The analysis was performed as described in Sahlmann 
et al.| ( |201 lb) , where a detailed description of the method can 
be found. Using the orbital parameters given by the RV analy- 
sis (Table |6]l, the IAD was fitted with a seven-parameter model, 
where the free parameters are the inclination /, the longitude of 
the ascending node Q., the parallax cr, and offsets to the co- 
ordinates (Aa*, A5) and proper motions (AyU„*, A/i^). A two- 
dimensional grid in / and Q. was searched for its global x^- 
minimum with a standard nonlinear minimisation procedure. 
The statistical significance of the derived astrometric orbit was 
determined with a permutation test employing 1000 pseudo or- 
bits. Uncertainties in the solution parameters were derived by 
Monte Carlo simulations. This method has proven to be reliable 
in detecting orbital signatures in the Hipparcos IAD ( Sahlmann 
et al.|201 1b a). Where the astrometric signal was not detected in 
the data, i.e. the derived significance is low, the Hipparcos obser- 
vations were used to set an upper limit to the companion mass 
foHowing the method described by Sahlmann et al. ( 2011b|a i. 

Table|5]lists the target identifiers and the basic parameters of 
the Hipparcos observations relevant for the astrometric analysis. 
The solution type indicates the astrometric model adopted by the 
new reduction. For the standard five-parameter solution it is '5', 
whereas it is '7' when the model included proper- motion deriva- 
tives of first order to obtain a reasonable fit, which is the case for 
HD167215. The parameter A^orb represents the number of orbital 
periods covered by the Hipparcos observation timespan and A^Hip 
is the number of astrometric measurements with a median preci- 



4. Results 

In this section we discuss the results for each target. We have 
classified the reported companions roughly according to their 
minimum mass and evidence for or against them being stellar 
companions. 

4.1. Planetary candidates 
HD22781 

HD22781 (HIP 17187) is a V=8.8 magnitude KO star located at 
d=33 parsecs from the Sun, with a stellar mass Ms - 0.75 Mq. 
The RV measurements obtained with SOPHIE (Fig. [T} reveal 
the presence of a sub-stellar companion with minimum mass 
McSin/ - 14 Mjup in a highly eccentric (e - 0.82), 528-day 
period orbit. Due to the effective temperature of the target, the 
velocities used to obtained the Keplerian fit are those obtained 
using the K5 mask. 

The residuals of a single-Keplerian fit exhibit a scatter of 
around 12 ms"', largely in excess of the typical measurement 
error, and too large to be explained by the effect of stellar jitter, 
estimated to be about 1 ms"' ((XogR'^j^ - -5.0). Furthermore, 
the residuals seem to exhibit a trend, with a maximum around 
BJD = 2455 100 and a subsequent decrease (see Fig.|2]). In view 
of this, the RV measurements were also fit using a Keplerian or- 
bit plus a quadratic drift. The best-fit parameters of this solu- 
tion are reported in Table |6] We note that the parameters corre- 
sponding to the Keplerian orbit are consistent at the 1-cr level 
with those derived from a single-Keplerian fit. Concerning the 
parabolic drift, we can obtain a rough estimate of the possi- 
ble body producing such a drift by assuming that its orbital pe- 
riod and the RV semi-amplitude correspond roughly to twice the 
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Table 6. Fitted orbits and planetary parameters. 



HD22781'' 



HD 156279 



HD 14651 



HD30246 



Orbital Parameters 



P [days] 

e 

K [m s-'] 

Tq (periastron) [B JD-2 454 000] 

Vo [kms-'] 
di [ms^Vyr] 
[ms^Vyr^] 
Fit Parameters 



528.07 ±0.14 
0.8191 ±0.0023 
-44.08 ± 0.56 
726.4 ±7.1 
881.40 ±0.12 
8.2508 ± 0.0042 
14.2761 ±0.0019 
-3.10954 



131.05 ±0.54 
0.708 ±0.018 
-95.8 ± 1.7 
578 ± 20 
1525.59 ±0.18 
-20.6432 ± 0.0050 



79.4179 ± 0.0021 
0.4751 ± 0.0010 
-151.92 ±0.21 
2595.0 ±3.6 
537.607 ± 0.045 

52.7244 ± 0.0030 



990.7 
0.838 ± 



t5.6 
0.081 



-65+"* 

^"-'•^-560 

1314.0!^i 
41.91!'* 



Typical photon noise 
Typical uncertainty" [m s"'] 
o-Q-c [ms-'] 
Derived Parameters 
/(Ms,Mp)[10-''Mo] 
a [AU] 

Min. astrometric signal [10^^ arcsec] 
Mc sin [Mj„p] 



M. 



c,up— lim 

M, [Me]" 



[Mo] 



2.61 
4.41 
8.36 

3.956 ± 0.066 
1.167 ±0.039 
22.62 ± 0.26 
13.65 ±0.97 
[0.23-0.31] 



2.83 
4.26 
9.08 

0.926 ± 0.034 
0.495 ± 0.017 
5.68 ±0.31 
9.71 ±0.66 
[0.35-0.41] 



3.17 
4.64 
5.26 

97.95 ± 0.41 
0.361 ± 0.013 
36.71 ±0.03 
47.0 ±3.4 
[0.28-0.4] 



4.35 
5.67 
14.09 



2.012 ±0.069 
48.8+1^^ 

[0.28-0:37] 



Orbital Parameters 
P [days] 

e 

ojn 

K [m s-'] 

To (periastron) [B JD-2 454 000] 
Vo [kms-'] 
di [ms^'/yr] 
d2 [ms-'/yr^] 
Fit Parameters 



HD92320 

145.402 ±0.013 
0.3226 ± 0.0014 
-22.99 ± 0.54 
2574.3 ± 6.6 
798.17 ±0.18 
-0.1418 ±0.0045 



HD137510 

801.30 ±0.45 
0.3985 ± 0.0073 
31.2± 1.1 
531.5 ±7.0 
187.8 ±2.0 
-6.4157 ±0.0081'' 



HD 16702 

72.8322 ± 0.0023 
0.1373 ±0.0017 
-102.44 ±0.56 
2423.7 ± 3.0 
983.45 ±0.12 
4.5602 ± 0.0022 



HD167215 



632+^60 
37+°-'^ 

"•-"-0.06 

140 ± 15 

1748+go 
-45 17+°-^'' 



Typical photon noise [m s ' ] 


3.31 


3.85 


3.32 


7.51 


Typical uncertainty" [m s"'] 


4.69 


5.32 


4.78 


8.16 


ctq-c [ms"'] 


9.07 


10.76'' 


7.85 


7.92 


Derived Parameters 










/(Ms, Mp) [10-" Mo] 


217.9 ± 1.7 


9.61 ±0.30 


104.41 ±0.39 


520+''''" 

-^^"-250 


a [AU] 


0.536 ±0.018 


1.880 ±0.064 


0.344 ± 0.012 


1 56+°'" 


Min. astrometric signal [10^^ arcsec] 


75.2 ±0.2 


86.3 ±0.9 


34.59 ± 0.06 


144+ [06 


Mc sin [Mj„p] 


59.4 ±4.1 


27.3 ± 1.9 


48.7 ± 3.4 


q9+S 


Mc.„p-ii,^ [Mo] 


[0.26-0.34] 


0.057 






[Mo] " 




[0.019-0.057] 


[0.35-0.45] 


[0.2-0.29] 



Notes. A second companion is probably present in the system (see 5 4. 1 1. The orbital parameters for HD22781 b might therefore change slightly 
once the two-Keplerian model is used. Including wavelength calibration error and 3 ms"' of systematic error. 99.7%-confidence interval. 
*'^* Systemic radial velocity of SOPHIE data. The values obtained for the other instruments are: 0.1587 ± 0.0082 km s"', 0.1352 ± 0.083 km s^'and 
0.045 ± 0.011, for McD, TLS and Lick data, respectively. '''' Dispersion for SOPHIE data. The values for other data sets are: 13 ms"', 20 ms~', 
and 16 m s^'for McD, TLS and Lick data, respectively. 



time between the zeros of the parabola, and to its maximum, 
respectively. In this way the fit coefficients imply a companion 
with minimum mass around 1.0 Mjup in a 9.2-year orbit. A two- 
Keplerian fit, on the other hand, produces a fit with half this 
period and a slightly smaller minimum mass (0.8 Mjup) mass 
for the second companion, and reduces the scatter in the resid- 
uals to about 7 ms"', which is comparable to the scatter level 
of the other targets. However, the fit semi-amplitude has a sig- 
nificance below 2-cr. We therefore decided to keep the solution 
with a quadratic drift, and note that this system warrants further 
observations in order to confirm the presence of a second orbit- 
ing body. We therefore warn the reader that the orbital elements 
reported in Table [6] should be taken as preliminary. 



The bisector velocity span of these observations is compati- 
ble with a constant (reduced smaller than unity) and a linear 
regression of the bisector as a function of RV yields a flat slope 
{b < 0.012, at 3-cr confidence level). These results hold for all 
spectroscopic masks employed, reinforcing thus the case for a 
sub-stellar companion. Additionally, no significant difference is 
observed in the semi-amplitudes obtained using different masks. 
The astrometric orbit derived from Hipparcos data has a low sig- 
nificance of 22 %, but we can set an upper limit of 0.31 Mq for 
the mass of the companion, which is increased by the low num- 
ber of 69 Hipparcos measurements with moderate precision and 
the high eccentricity of the system. If the secondary is a slow- 
rotating star, this constraint can be slightly improved by means 
of the CCF simulations described in § 3.5. 1 which give an upper 
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Fig. 1. Radial velocity measurements obtained with SOPHIE for all objects studied in this paper. In the left column, the measure- 
ments are plot as a function of time, while on the light column, they have been folded to the orbital phase. The solid curve represents 
the best fit to the data and the lower panels on the left column show the residuals to the fit. 



limit of 0.23 Mq. In the case of a fast rotating star, the upper 
limit is similar to that from astrometry. 

We conclude that the companion to HD22781 is most likely 
a sub-stellar object, which could have formed either by core- 
accretion (Mordasini et al. |2009[l or by gr avitational collapse. 
The term "superplanet" used by |Udry et al!] ( |2002[ ) for objects of 



similiar mass seems adequate (see The detection of a poten- 
tial second companion -albeit with low significance- also points 
to the planetary nature of HD22781 b, although it is known that 
circumbinary planets exist ( |Doyle et al.||201 1} . The dynamical 
stability of this scenario could be promptly tested by numeri- 
cal simulations, but this is outside the scope of this paper. In 
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Fig. 1. Cont. For HD137510, the red, green and blue points correspond to data from McDonald, TLS, and Lick, respectively. 



any case, this system warrants additional RV observations over 
a longer time period. 

HD1 56279 

HD156279 (HIP 84171) is a G6 dwarf star with a mass of 0.93 
Mq. It is located 35 pc away from the Sun and has a V-band 
magnitude of 8.08. The RV variations observed with SOPHIE 



imply the presence of a planetary-mass companion (Mc sin / = 
9.71 ± 0.66 Mjup) in a 130-day eccentric orbit. The best-fit orbit 
is shown in Fig.jTjand the parameters are reported in Table|6] 

Bisector velocity span measurements for this star do not ex- 
hibit any significant variation in time, nor any dependence with 
RV, and present a scatter of 6.5 ms '. The CCF simulations 
(Fig. [3]l put an upper limit of roughly 0.35 Mq, although so- 
lutions with higher mass are possible for rotational velocities 
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HD156279 



HD30246 




9 10 11 12 

FWHM2 [km s '] 

Fig. 3. Confidence regions for the mass of the secondary and 
the FWHM of its CCF, issued from numerical simulations of 
HD 156279. The dashed line is the upper limit derived from the 
astrometric analysis of Hipparcos data. At the 3-cr level, the 
mass of the secondary is constrained by these simulations to 
< 0.35 Mq. Astrometric data discard the high-mass family with 
FWHM ~ 7 km s"' permitted by the bisector data. 



around the resolution Umit of SOPHIE (FWHM » 6-7 kms 
depending on spectral type). However, this family of solutions 
is discarded by Hipparcos measurements, for which the derived 
upper mass is 0.41 Mq. Additionally, no mask effect is detected 
at a significant level, which also supports the hypothesis of a 
sub-stellar object orbiting HD 156279. 

HD156279 b is the least massive object reported in this pa- 
per, and may have been formed as a planet. Its minimum mass 
is below the smallest deuterium-burning mass determinations 
HD 156279 b is therefore established as 



(Spiegel et al. 2011 



a planetary candidate. 



4.2. Brown dwarf candidates 
HD14651 

HD14651 (HIP 11028) is a V=8.3 magnitude star, with a mass 
Ms = 0.955 Mq, measured as described in Sect. |3.1| The 
Keplerian fit to the RV measurements obtained with SOPHIE 
(Fig. [T} indicates that this star possesses a companion of mini- 
mum mass Mc sin / - 47 Mj in a 79.4-day orbit, with an eccen- 
tricity e - 0.48. 

Evidence that the companion is indeed a sub-stellar object 
comes from the fact that the amplitude of the orbital fit does not 
depend on the numerical mask used to compute the CCF, and 
from the bisector analysis, which is compatible with a constant 
(reduced - 0.6), and does not exhibit any correlation with the 
RV measurements. Indeed, the Pearson and Spearman correla- 
tion coefficients between the RV and the bisector velocity span 
are -0. 1 1 and -0. 14, respectively, and the slope obtained by linear 
regression is consistent with zero at the 1-cr level. Additionally, 
we do not detect an orbital signature in the Hipparcos astrom- 
etry, the significance being 17 %. We can set an upper limit of 
0.72 M0 for the mass of the companion, which is high partly 
due to the low number of Hipparcos measurements (A^Hip - 59) 
with moderate precision and the considerable eccentricity of the 
system. On the other hand, the CCF simulations constrain the 
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Fig. 5. Confidence regions for the mass of the secondary and 
the FWHM of its CCF, issued from numerical simulations of 
HD30246. At the 3-cr level, in addition to the upper limit at 
around 0.33 Mg, there exists a solution with ~ 0.45 Mg. 



mass of the secondary to be smaller than 0.4 Mg if the star is 
rotating with projected velocity v sin / a: 10 km s"' , and smaller 
than 0.28 Mg if V sin/ < 3kms-'. 

We conclude that the companion to HD 14651 is a bona 
fide brown dwarf candidate, that requires further observations to 
pinpoint its true nature. For example, astrometric observations 
from the future GAIA mission, or near-infrared spectroscopy 
that could detect the peak of the companion, and therefore mea- 
sure the mass ratio between primary and secondary (see 

HD30246 

HD30246 (HIP 22203) is a V=8.3 magnitude Gl star located at 
5 1 parsecs from the Sun, with a measured mass - 1 .047 Mg. 
The RV measurements obtained with SOPHIE exhibit a peak-to- 
peak variation of 2 km s"' , that can be explained by the presence 
of a Mc sin / - 56 Mj companion in a 990-day orbit, with an 
eccentricity e = 0.85. HD30246 belongs to the Hyades cluster 
( |Griffin et al.| |1988 ) and has been the object of several studies, 
discussed in part below. It is also the most active star in the sam- 



ple, compatible with its probably young age (Holmberg et al. 
|2009l ). 

For this target, observations do not cover one of the extrema 
of the RV curve (see Fig.[T]i. As a consequence, the Monte Carlo 
distribution of the minimum mass of the companion has a long 
high-mass tail, and a slightly bimodal behaviour, allowing for 
a family of solutions corresponding to a low-mass star (Fig. [4]i. 
The distributions of the semi-amplitude K and the orbital eccen- 
tricity e exhibit a similar feature, with the high-mass solutions 
coiTesponding to high-eccentricity orbits. Since the value of the 

statistics is not systematically higher for any value of the ec- 
centricity or the semi-amplitude K, the identification of the cor- 
rect solution is hindered. We note, however, that less than 14% 
out of the 5000 Monte Carlo realisations of the data produce an 
orbital solution with minimum mass above 80 Mjup. The best fit 
orbital solution is shown in Fig. [T] and the obtained parameters 
obtained using mask G2 are reported in table |6] 

The bisector velocity span is compatible with a constant, ex- 
hibiting a dispersion of around 16 ms '. Also, no relation is de- 
tected with the RV above the 1.8-cr level {b = -0.0059 + 0.0034), 
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Fig. 4. Left: Minimum mass distribution for HD30246, obtained from 5000 Monte Carlo simulations. The bimodal behaviour is 
clearly seen, and allows for stellar-mass solutions. Right: CoiTelation between orbital eccentricity and RV semi-amplitude from the 
same Monte Carlo simulations, showing that high-mass solutions correspond to very eccentric orbits. 



and the amplitudes of the RV variation measured with differ- 
ent masks are in agreement {Kq2 - 1995.1 ± 4.1 ms"', Kf^^ - 
1996.8 ± 4.6 ms-i, Km4 = 2005 + 16 ms"'). Additionally the 
large scatter in the fit residuals (around 13 ms"') canbe partially 
explained by the stellar jitter of the host star, which is estimated 
to be around 10 ms~ ', given the median value of activity index 
(log/?^^) = -4.50 (|Santos et al.]]2000|l. Also, the analysis of 
the Hipparcos astrometfy proviTie orbits with low significance 
(87.6 %, or 1 .5 cr), meaning that the orbital signature is not de- 
tected. However, because the orbit is not covered and in addition 
is highly eccentric (e ~ 0.8), we refrain from deriving an up- 
per limit to the companion mas^ On the other hand, the CCF 
simulations described in 5 3.5. 1| allow us to set an upper limit 
at the 3-cr level of 0.28 Mq (respectively, 0.37 Mq) for simula- 
tions in the slow-rotating (respectively, fast-rotating) extreme of 
the explored parameter range. Additionally, a possible solution 
is found for Mc in the interval [0.42 - 0.51] Mq, with rotational 
velocity v sin / between 3 km s"' and 4.5 km s"' , at the 3-cr level 
(see Fig.[5]l. 

HD30246 has been reported to be a spectroscopic binary 
( [Griffin et al.| 1988'), based on 13 RV measurements with a preci- 
sion of around 700 m s ' . Additionally, signs of the spectrum of 
the secondary were found in 4 infrared spectra obtained with the 
NASA Infrared Telescope Facihty by [Bender & Simon] ( [2008l l. 
However, according to these authors, the detection of a sec- 
ondary peak in the CCFs of this target is challenging, and leads 
to poorly-constrained velocities for the secondary (see §5.3 and 
Figure 4 of their paper). The derived mass ratio q - 0.06 ± 0.04 
implies a mass Mp - (66 ±45) Mjup for the companion, when the 
stellar mass derived in this paper is considered, - 1 .05 + 0. 1 1 
This is in agreement with our determination for the mini- 
mum mass of the companion. Additionally, the orbital parame- 
ters of HD30246 b, obtained by Stefanik & Latham and reported 
by Bender & Simon ( 2008| l, agree with the parameters we de- 
rived from SOPHIE data. Furthermore, K-band speckle interfer- 
ometry ( [Patience et al.[ 1998 ) provides a limiting magnitude dif- 
ference for companions at 0.05" from HD30246 of AK - 4.6. 
Considering the V - K relations for main-sequence stars (Cox 



Even when taking the obtained solutions at face value, the derived 
mass hmits of ~50 - 400 Mj would not allow us to establish the nature 
of the companion. 
^ A 10% error in the stellar mass is used, rather than the formal error 



2000[) and the mass-luminosity relations from Delfosse et al. 
( 2000| l, this translates roughly to a maximum mass of 0.19 Mg 
for companions to HD30246 at this orbital distance, which is 
slightly larger than the distance obtained from the Keplerian fit 
to the data, a - 0.04". 

We conclude that the initial claim of binarity by|Griffin et al. 
( 1988| l warrants revision, and that the remaining facts support a 
sub-stellar candidate companion to HD30246. The upper limit 
for the mass is however in the stellar regime, and therefore fur- 
ther astrometric observations will be needed to decide unam- 
biguously on the nature of this object, as well as a continuing 
RV monitoring, in order to better constrain the orbital parame- 
ters. 

HD92320 

HD92320 (HIP 52278) is a V=8.38 magnitude star located at 41 
parsecs from the Sun. The 16 RV measurements obtained since 
late 2006 are fitted by a 145 -day Keplerian orbit with a semi- 
amplitude of about 2.6 kms and an eccentricity of 0.32. The 
spectral analysis coupled to the comparison with evolutionary 
models yields a stellar mass of - 0.916 Mq, and a companion 
minimum mass Mc sin / = 59.5 Mjup. 

The Keplerian fit residuals exhibit a weak correlation with 
pupil illumination proxy d^v (see 5 3.4.1i, with Pearson and 



Spearman correlation coefficients -0.60 and -0.66, respectively. 
The trend is given mainly by the four points with largest absolute 
residuals (see Fig. [6]). These measurements do not exhibit un- 
usual signal-to-noise ratios, CCF contrast or FWHM. We have 
therefore concluded that the relation is real and consequently, 
the detrended measurements were used to obtain the orbital fit 
plotted in Fig. [T] and the parameters reported in table |6] 

The lack of cotTelation between bisector velocity span and 
R\0 together with the absence of mask effect argues for a 
substellar companion as the source of the periodic RV varia- 
tions. Additionally, we did not detect an orbital signature in the 
Hipparcos data, but we could set an upper limit of 0.34 Mq for 
the mass of the companion. The CCF simulations described in 



3.5.1 provide a similar upper limit to the mass: at the 99.7% 



confidence level, the upper limit to the true mass of the com- 
panion ranges from 0.26 Mq, for slow rotators below the detec- 



reported in Table|4] 



' The Pearson's and Spearman's coefficients are -0.083 and -0.144, 
respectively 
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HD92320 ^ ^ HD137510 HD16702 




Fig. 6. Residuals of the Keplerian orbital fit to the uncorrected data as a function of pupil illumination indicator 6rv for the three stars 
in our sample showing a significant (r > 0.5) correlation. The line is a fit to the data, and was used to correct the radial velocities, 
which leads to a reduced scatter. 



tion limit up to 0.35 M© for the fastest rotating stars simulated 
(vsin/ = 10 kms '). The companion to HD92320 is therefore 
established as a brown dwarf candidate that merits further obser- 
vations. GAIA measurements should easily detect the reflex mo- 
tion induced by this object (^, and near-infrared spectroscopy 
might permit detecting the spectrum of the secondary. 



HD137510 



HD137510 (HIP 75535) is a high-metaflicity 1.36-Mo sub-giant 
late-F star, located 42 pc away from the Sun. Its V-band mag- 
nitude is 6.26, making it the brightest star in our sample. A [Quirrenbach 1 20 1 if 
brown dwarf companion candidate with a minimum mass of 26 
Mjup was discovered by jEndl et al.| (j2004), based on observa- 
tions acquired at the McDonald Observatory (McD) in Texas, 
USA and at Thiiiinger Landessternwaite Tautenburg (TLS), in 
Germany. Later, [Butler et al. ( 2006 1 published an independent 
orbital solution based on 10 RV measurements obtained at Lick 
Observatory, which led to a slightly less massive object of 22.7 



The scatter in the bisector velocity span measurements is be- 
low 16 ms', and produced mainly by one point, which if left out 
reduces the scatter to around 12 ms"', consistent with the typi- 
cal error of the RV data (~5 ms '). Additionally, no correlation 
with the RV measurements above the 1.5-cr level is present, and 
there is no significant difference in the RV amplitudes measured 
through different numerical masks. The results of the CCF sim- 
ulations are compatible with a substellar object at the 3-cr level, 
yet another solution with Mc ~ 0.5 Mq is also possible. This 
last solution is discarded by the astrometric analysis. Indeed, the 
substellar nature of HD 1375 10 b has already been established 
byjSahlmann et al. (201 lb I and later confirmed by [Reffert & 



Sahknann et al. ( 2011b| l derived an upper 



Mjup, in a mildly less eccentric orbit. Thirteen new measure- 
ments of this star were obtained with SOPHIE, which permit 
improving the system parameters. Indeed, the determination us- 
ing the combined datasets has reduced uncertainties, particularly 
in the orbital period and in the eccentricity. 

The SOPHIE observations were corrected as described in 
^3.4.1 1 for the seeing effect. The Pearson's and Spearman's co- 
efficients between the residuals and (5rv are -0.67 and -0.70, 
respectively (Fig. |6|l. The dispersion of the residuals of the 
Keplerian fit to the SOPHIE data is improved by this procedure: 
from around 8.6 m s"' for the uncorrected data to 4.6 m s"' . The 
fit to the combined data of all four instruments is also improved, 
albeit less drastically: the scatter of the residuals is reduced from 
15.2 m s^'to 14.4 m s ' . In Table[6]we report the parameters ob- 
tained from the combined fit of all four datasets, and the best fit 
curve is shown in Fig.[T] 

Concerning the stellar parameters, we find an effective tem- 
perature Tgff = 6130 ± 50 K, in good agreement with the value 
by Lebr eetaL] ( [T999] l, and about 160 K (3 .1-cr)and 230 K (2.8- 
cr) larger than the value reported by Valenti & Fischer ( 2005 ) 
and j Endl et al. (2004i, respectively. These discrepancies might 



companion-mass limit of 64.4 M jup using H ipparcos astrometry 
and the spectroscopic elements of [Butler eFal. (2006). On the ba- 
sis of the updated RV solution presented here, we find a similar 
result: the astrometric signature is not detected, but because of 
the updated orbital elements and the slightly lower stellar mass, 
we obtain an updated upper limit of 59.5 Mjup for the mass of the 
companion, which confirms that HD 1375 10 b is a brown dwarf 
companion orbiting an F-star. 



be due to unaccounted systematic errors in the different methods 
and models used. The metallicity found as described in Sect. |3.1 
([Fe/H] - 0.38 + 0.13) is in agreement with previous determina- 
tions, and agrees also with the metalhcity value obtained from 
the CCF as described in §3.1| 



4.3. Stellar companions 
HD16702 

HD16702 (HIP 12685) is a V=8.3 magnitude GO star located at 
45 parsecs from the Sun. The orbital fit to the RV data (Fig. [TJ 
indicates that this star possesses a companion in a 72.8-day orbit, 
with an eccentricity e = 0.14 and minimum mass sin / - 48.7 
Mjup, given a stellar- mass Ms = 0.98 Mq. 

The residuals of the keplerian fit exhibit a clear coiTelation 
with Srv (see Figure. |6]l, which allowed us to correct the data, as 
described in ^3.4.1| When this is done, the scatter in the fit resid- 
uals is reduced from 12 ms"' to the final value of » 8 m s"' . The 
final fit parameters are listed in Table |6] and the best-fit curve is 
plotted in Fig.[T] 

The bisector analysis for this target shows a significant corre- 
lation between bisector velocity span and radial velocity, shown 
in Fig.|7] A linear regression to the data shows that the slope b is 
significantly different from zero, b = -5.1 + 1.3, i.e. a detection 
at the 4-cr level. Additionally, the amplitude of the RV variations 
depends on the spectral mask used to compute the CCF. Indeed, 
we find Kc2 = 2423.7 + 3.0 ms"', Kfcs = 2412.0 + 5.2 ms"', 
and Km4 - 2344.0 ± 9.5 ms"', which implies a 8-sigma dif- 
ference between mask G2 and mask M4. These facts indicate 
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- HD 16702 9 



RV [km s '] 

Fig. 7. Bissector velocity span and radial velocities for HD 16702 
(corrected for the seeing effect). The solid line represents the best 
fit to the data and the dashed lines indicate the 3-cr confidence 
region. A correlation is detected to 4-cr significance, which hints 
to the presence of an additional blended star in the fiber of the 
spectrograph. 



that the companion around HD 16702 might be a star in a low 
inclination orbit. The spectra obtained would then consist of the 
addition of the spectrum of the solar-like primary star and the 
secondary star, probably a low-mass star, moving in anti-phase 
in velocity space. When correlated with the M4 mask, the rela- 
tive weight given to the spectrum of the secondary is increased, 
and the resulting CCF is then skewed towards the velocity of the 
secondary, which produces a decrease in the measured amplitude 
of the RV variations. 

The astrometric analysis and CCF simulations confirm this 
scenario. We detect the astrometric signal in the Hipparcos data 
with a significance of 2.5 cr (98.7 %) and measure a semimajor 
axis of 2.9 + 0.6 mas. We see the system almost face-on with 
an orbital inclination of i — 6.6;|;[ j degrees, corresponding to 
a companion mass of Mc = 0.55 + 0.14 M0 and we measure 
the longitude of the ascending node of Q = 146^^*^ °. The co- 
ordinate and proper motion offsets are Aa* = -0.4 + 0.7 mas, 
A6 = 0.0 + 0.5 mas, A/i„* = -1.8+ 1.2 mas yr"', Afis = 0.6 ±0.8 
mas yr and the updated parallax ism = 20.9 + 0.8 mas, which 
is compatible with the original value. The sky-projected orbit of 
HD16702 and the Hipparcos measurements are shown in Fig. [8] 
The values of 2.5 cr for the significance and of 6.5 for the astro- 
metric S/N for this orbit indicate that we are close to the limit 
of the detection capabilities of Hipparcos. However, our astro- 
metric analysis rales out the possibility of the companion being 
a sub-stellar object. Indeed, at the 3-cr level, the lower and upper 
companion mass limits are 0.19 Mg and 1.05 Mq, respectively. 

The CCF simulations permit further constraining the com- 
panion mass. In Fig.|9]we show the c onfiden ce regions obtained 
from the;\f^ surface, as described in §3.5.1 We deduce that the 
maximum allowed value, at 99.7 % confidence level, range from 
0.35 Mq to 0.45 Mq, depending on the FWHM of the secondary. 
Interestingly, due to the significant correlation between bisec- 
tor and radial velocity shown in Fig. [7] also a lower limit of 
about 0.25 Mg can be put to the mass of the secondary. On the 
other hand, no strong constraints can be put on the rotational 
velocity of the secondary. We note that the best-fit point has 

= 14.4, comparable with the number of observations (N=22), 
which means that the confidence regions are a rigorous way to 
estimate the mass of the secondary. 
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Fig. 8. Top: Astrometric orbit of HD 16702 projected on the sky. 
North is up and east is left. The solid red line shows the model or- 
bit, which is orientated counter-clockwise, and open circles mark 
the individual Hipparcos measurements used for the astrometric 
fit. Blue filled circles show normal points for each satellite orbit 
number and are used for display only. Bottom: O-C residuals for 
the normal points of the orbital solution (filled blue circles) and 
the five-parameter model without companion (crosses). 



HD 16702 is therefore a binary star, with a mass ratio ranging 
from q « 0.25 to q ^ [0.35 - 0.46]. See ^for further discussion 
on this object. 



HD167215 

The spectral analysis of HD 1672 15 (HIP 89270) shows that 
this object is a low-metallicity 1.15-Mq late-F star. It is located 
79 pc away the Sun, and the RV measurements obtained with 
SOPHIE, although inconclusive due to their incomplete phase 
coverage, point at a 600-day period variation with an amplitude 
of 2 kms"'. 

The Monte Carlo distribution of minimum masses (see 
Fig. 10 1 exhibits a clear bimodal behaviour, with a lot of power 



on the stellar-mass range. Indeed, only about 25% of the Monte 
Carlo realisations of the data are fitted to a minimum mass below 
80 Mjup. Therefore, the orbiting object is most probably a low- 
mass star. Just as for HD30246, the high-mass solutions corre- 
spond to high eccentricity orbits. Also, the period is well corre- 
lated with the semi-amplitude K (Fig. 



lOi 



The bisector velocity span is presented as a function of RV in 
Fig. 1 1 A constant fit to the data gives;^^^ - 23, for 13 degrees of 



13 



Diaz et al.: Massive companions in the planet-brown dwarf boundary. 




40 60 ao 100 120 140 160 180 200 

Minimum mass [Mjupl 




RV semi-amplitude [km s ^] 



Fig. 10. Left: Minimum mass distribution for HD167215, obtained from 5000 Monte Carlo simulations. Only about 25% of the 
solutions lie below the 80-Mjup limit. Right: Correlations between the RV variation amplitude and the orbital eccentricity (upper 
panel) and the orbital period (lower panel). 
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Fig. 9. Confidence regions for the mass of the secondary and 
the FWHM of its CCF, issued from numerical simulations of 
HD16702. At the 3-cr level, a strong constraint can be put to the 
mass of the secondary. 



freedom, and a linear regression gives - 21.6, for 12 degrees 
of freedom. An F-test and a run test (see [Frodesen et al.|[T979l 
§14.6.2) fail at rejecting the null hypothesis at 5% significance 
level. Additionally, no evidence of mask effect is detected, which 
means that no significant amount of light from the putative sec- 
ondary star enters the spectrograph. Finally, we did not detect 
an orbital signature in the Hipparcos data, since the derived or- 
bit has a significance of only 56 %. We were nevertheless able 
to set an upper limit of 0.29 M0 for the mass of the compan- 
ion. When the RV and astrometric analyses are combined with 
the CCF simulations (Fig. 12 1, the mass of the secondary is con- 
fined to the wedge-like region of Fig. 12 delimited by the upper 
limit set by astrometry (0.29 Mq, shown as a dashed line) and 
the 99.7 % contour level below it. The secondary is likely to 
have a projected rotational velocity below vsin/ ^ 5.3 kms"', 
which is the value coiTespoding to FWHM » 10 kms The;\f^ 
of the best fit is 1 1.3, for 14 observations. Regardless of the low 
significance of the bisector velocity span variation, the CCF sim- 
ulations place the mass of the secondary in the stellar domaine. 
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Fig. 11. Bisector velocity span as a function of RV measure- 
ments for HD167215. The curve is the parabolic fit to the data. 



indicating that the companion to HD 1672 15 is likely a low-mass 
star, with mass ratio between ^ = 0.17 and 0.25. 



5. Discussion 

5.1. The SOPHIE objects 

The objects reported in this paper span a wide range of min- 
imum masses, from 9.5 Mj^p (HD 156279 b) up to 92 Mj^p 
(HD167215B), and are therefore useful to perform a compar- 
ative study of the properties of massive planets, brown dwarfs 
and low-mass stars. In particular, the object orbiting HD137510 
has been established to be a brown dwarf ( jEndl et al.||2004[ 
I Sahlmann et al.|[2011b[ and this work), and the companions 
to HD 16702 and HD 1672 15 have been identified as low-mass 
stars. Out of the remaining 5 objects reported here, two have 
minimum masses below or close to the deuterium-burning limit, 
and we have labelled them as planet candidates (HD22781 b 
and HD 156279 b), while the other three are brown dwarf can- 
didates. Except for HD92320 b, which has a minimum mass 
slightly larger, these candidates fall in what has been identified 



by Grether & Lineweaver ( 2006 ) as the driest part of the brown 



43^23 Mjup. This implies that these can- 



dwarf desert, at Mc 

didates are rare and therefore worthy of further observations in 
order to pinpoint their true mass. In Table |6] we report the min- 
imum astrometric signal expected for the detected candidates. 



14 



Diaz et al.: Massive companions in the planet-brown dwarf boundary. 




60 



o RV candidates 

■ SOPHIE 

• RV + Astrometry 

• Transit 



Orbital Period [d] 



u 0.6 



o RV candidates 

■ SOPHIE 

• RV + Astrometry 

• Transit 




Orbital Period [d] 



Fig. 13. Upper: Mass-period diagram for all objects with (minimum) masses between 9 and 90 Mjup, orbiting solar-type stars with a 
period shorter than 1 x lO'* days. Lower: Eccentricity-period diagram for same sample of objects. In both panels, the objects reported 
in this paper are shown as red squares. The other coloured symbols represent companions for which the true mass is known, either 
by astrometry (blue circles) or because they transit their host star (green circles). The rest of the objects are candidates detected by 
RV surveys (white circles), and the plotted mass should be interpreted as a lower limit. Due to the large range of objects represented, 
for some of them, the error bars are smaller than the symbol used in the plot. 



defined as the semi-amplitude of the periodic motion of the star 
projected on the plane of the sky, assuming / - 90 degrees. The 
future space mission GAIA, having an astrometric precision be - 



low 1.0x10 arcsec for blight stars such as these (Prusti 



2011 



should easily detect the orbital signature of all three of these can- 
didates. 

HD 16702 has been identified as a stellar binary system, with 
a period of 72.8 days and a mass ratio q lying in the interval 
[0.25-0.46] at the 99.7% confidence level. Considering the mass 
of the companion to be - 0.35 Mq, the magnitude differ- 
ences in the V and K bands imply a contrast of 3 x 10"^ and 
5.6 X 10"^, respectively, using the apparent ^T-band magnitude 
from the 2MASS catalogue (Cutri et al. 2003) . The measured 
projected separation for this system is 8.6 mas. It will not be 
easily resolved with the current adaptive optics facilities. On the 



other hand, interferometric observations should be able to re- 
solve the system. 

For HD 1672 15, we find a binary system at a period of over 
600 days, and a mass ratio in the interval [0.17-0.25]. If we con- 
sider that the secondary has a mass of 0.25 Mq, inversion of the 
mass-luminosity relations by |Delfosse et al.| ( |2000[ ) provide the 
absolute magnitudes My - 12.41 and Mf^ - 7.28, in the V and 
/Tband. The primary star has apparent magnitudes my = 8.1 and 
niK - 6.69 (Cutri et al.|2003i ), which give contrasts of 2.8 x 10 
and 8.6 x 10"^ in the V and K band, respectively. The measured 
mass for the companion implies that the inclination angle / is 
roughly 20 degrees, and the projected separation is 26.5 mas. 
This system will also be difficult to resolve with current adaptive 
optics instruments, but it might be a suitable target for interfero- 
metric observations. 
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Fig. 12. Confidence regions for the mass of the secondary and 
the FWHM of its CCF, issued from numerical simulations of 
HD167215. When combined with the astrometric analysis, the 
true mass of the companion is constrained to the region below 
the astrometric upper limit (dashed line) and the 99.7% contour 
below it. 



The brown dwarf orbiting HDl 375 10 has been discovered by 
Endl et al. ( 2004) and confir med later by Sahlmann et al. ( 2011b| l 
and Reffert & Quirrenbach (201 1 1. The minimum mass reported 
in this paper based on the analysis of all available datasets is 
compatible with previous determinations, and the analysis of the 
SOPHIE data alone also produces a compatible result. If we con- 
sider that the true mass of the brown dwarf is 30.9 Mjup, i.e. the 
expected value assuming random orientation of the system, and 
including the limit on the inclination provided by the astrometric 
analysis (sin ; > 0.4588), then linear interpolation of the model 
grids by Baraffe et al. (2003 ) give eff'ective temperatures between 
800 K and 870 K for the two limiting values obtained for the age 
of this system, 3.0 Gyr and 2.3 Gyr. This identifies this object 
as a T-type dwarf, which should present prominent methane and 
water features in its near-infrared spectrum (see, for example, 
|Kirkpatrick|[20 05). However, being about 21 mas from its host 
star, which outshines it by over 15 magnitudes in the H band, 
obtaining an spectrum of this object will prove challenging. 

HD22781 exhibits excess variation in the residuals of the 
single-Keplerian fit, that might indicate the presence of a sec- 
ond companion in the system. A tentative estimation of the mass 
and period of this object results in a minimum mass of about 1 
Mjup in a 9 years period orbit. Although this would be an un- 
usual system, with the more massive planet in the inner orbit 
and a lighter companion farther away, this configuration is by no 
means completely unheard of. For example, HD20226 ( Correia 
et al.||20()5 Udry et aLl|2002 1, present a similar configuration, 
with an inner massive planetary candidate in a 255-day orbit, and 
an outer companion close to the 5:1 resonance. Other examples 
of systems of this kind are the three-planet system H1P14810 
( [Wright et al. 2009 2 0071), HP 1 34987 Pones et al.|[20T0l i, and 
] ID200964 Pohnson et al.|201 l| l. 



5.2. Formation and evolution of sub-stellar companions 

In this section, we put in context the results obtained by consid- 
ering similar objects reported in the literature. We do not attempt 
to perform here a complete and thorough statistical analysis of 



the available data, which is defer to a future work, but only to 
present a few relevant characteristics of sub-stellar massive com- 
panions that might shed some light on mechanisms involved in 
their formation and evolution. 



5.2.1. Orbital period 

The SOPHIE objects cover a range of orbital periods from little 
over 70 days for HDl 6702 b up to 990 days for HD30246 b, 
or roughly 2 AU. The position of these objects in the mass- 
peiiod diagram is shown in the upper panel of Fig.^3l to gether 
with other similar objects presented in the li teratur^ S a hlmann | 
|et al.|201 IbHSozzetti & Desidera|20T0l and references therein). 
HD16702B and HD167215B are not included since their true 
mass is above the imposed limit of 90 Mjup. Note that none of 
our objects have been detected in short-period orbits. Indeed, it 
seems evident in Fig. 13 that the number of objects with masses 
above ~10 Mjup increases with orbital period, even though RV 



and transit observations are biassed towards discovering objects 
in short-period orbits. This might be explained by considering 
that these objects are formed far from the central star by grav- 
itational collapse, reaching the mass needed to open a gap in 
the disk ( Ward|1997) l before start migrating inwards by interac- 
tions with the disk. Since then migration will be of type 11 (Ward] 
[1997), which has longer timescales, these objects might not have 
time to reach closer orbits before the disk dissipates. 



5.2.2. Orbital eccentricity 



In the lower panel of Fig. 13 we plot the orbital eccentricity of 



the SOPHIE candidates as a function of orbital period. It can 
be seen that all the reported objects exhibit a considerable or- 
bital eccentricity. Ribas & Miralda-Escude ( 2007| ) report an in- 
crease in the orbital eccentricity of massive planetary objects, 
when compared with less massive objects. These authors argue 
that this might be indicative of different formation mechanisms, 
and propose that while low-mass objects form by core accre- 
tion in protoplanetary discs, more massive planets might form by 
gravitational collapse, either of a pre-stellar cloud or of a grav- 
itationally unstable circumstellar disk. Since high-mass planets 
would then need to form far from their stars, the fact that those 
observed at shorter distances (typically a few AU at most) have 
higher eccentricities than low -mass planets is, according to these 
authors, due to the large radial migration they have undergone. 
Indeed, numerical simulations ( Papaloizou et al.||2001 Nelson 
etal.|2000 ) have shown that eccentricity is excited on massive (> 
20 Mjup) migrating objects by means of interaction with the cir- 
cumstellar disc, but not so for less massive objects. On the other 
hand, interactions at disk resonances can damp eccentricity (see 
Tremaine & Zakamska|2004l and references therein). Therefore, 
although current data seem to agree with the idea that massive 
companions have larger eccentricities, it is not clear whether this 
can be explained by invoking a larger radial migration. 

5.2.3. Metallicity 

Planetary companions with masses comparable to that of Jupiter 
are known to occur more frequently around metal-rich stars 
( [Santos et al.|200l[|Fischer & Valenti|2 005), which is considered 
strong support for the core-accretion formation scenario ( Pollack] 



' ^ Objects witli m inimum mass above 9 Mjup and below 100 Mjup 
orbiting solar-type s tars (O.V Mp, < M < l.tiM^), with orbital periods 
shorter than 10"* days. 
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let al.|1996||Alibert et al.|2005| l. On the other hand, stell ar com- 
panions are more common around metal-poor primaries (Grether 
& Lineweaver 2007; Raghavan et al. 2010), a trend that is more 
evident for stars with (B - V)1>1)J625 . Raghavan et al. (2010") in- 
terpret this (tentative) result by suggesting that lower-metallicity 
clouds might be more likely to fragment to form binary stars, as 
is revealed in numerical simulations. If this is so, and if brown 
dwarf companions form as stars, then we might expect them to 
show a similar trend. 



Another trend reported by Ribas & Miralda-Escude ( 2007 1 
is that the metallicity of the host star decreases with companion 
mass. This is interpreted as a hint that two different mechanisms 
act in the formation of extrasolar planets, as discussed above, 
and that massive objects form as easily around metal-rich stars 
than around metal-poor ones. Of course, this statement assumes 
that the sample of target stars out of which extrasolar planets are 
detected is uniform in metallicity. However, since massive com- 
panions are rare and RV surveys usually have biased samples, 
the trend observed by Ribas & Miralda-Escude ( 2007 1 might not 
be real. The actual relevant measurable to study this issue is the 
occurrence rate of detected objects as a function of the metal- 
licity of their host stars. However, given the small number of 
detected brown dwarfs around solar-type stars, assessing the re- 
ality of any such trend is challenging. 

6. Summary and Conclusions 

We have reported the discovery of 7 new sub-stellar candidates 
orbiting solar-type stars with minimum masses above ~ 10 Mjyp. 
Out of these, we have managed to identify two objects as low- 
mass stellar companions, based on astrometric measurements 
and numerical simulations of their cross-correlation functions. 
Additionally, the orbital parameters of a known brown dwarf 
companion, namely HD137510 b, have been improved upon us- 
ing radial velocity measurements obtained with SOPHIE. The 
three confirmed objects are too faint and too close to their host 
stars to be easily detected by current adaptive optics techniques, 
since they lay within the diffraction limit of the biggest avail- 
able telescopes. On the other hand, the brown dwarf candidates 
should be promptly confirmed or rejected with GAlA's micro- 
arc second astrometry. 

Different techniques were employed with the aim of discern- 
ing low-mass stars in low-inclination orbits from actual sub- 
stellar objects. To do this, we employed not only traditional RV 
diagnostics, such as the mask effect and the bisector analysis, but 
we also analysed astrometric data from the Hipparcos mission, 
and developed a new technique to fully exploit the information 
present in the bisector velocity span. This latter technique has 
allowed us to set upper limits to the masses of all the candi- 
dates detected here, even surpassing in some cases the possibili- 
ties of current astrometric data, and has allowed the confirmation 
that one of the candidates was in fact a low-mass star. In some 
cases, the CCF simulations and the Hipparcos astrometry turn 
out to be complementary, as for HD167215. The effectiveness 
of each technique depends strongly on the orbital parameters, 
the number of measurements and the phase coverage, and we 
have therefore refrained from studying this here (see Santerne et 
al., in preparation). Improvements of the CCF simulations used 
in this method are currently underway, and should provide fur- 
ther possibilities for this promising technique, such as validating 
planetary candidates discovered by transiting surveys that show 
no variation in their RV or bisector velocity span measurements, 
as is the case for the recently-announced CoRoT-22 b (Moutou 
et al., submitted). 



Concerning the question of the dividing line between planets 
and brown dwarfs, if a mass-dependent classification is adopted, 
motivated by the observed mass distribution, such as the one pro- 
posed by Sahlma nn et al.| ( |2011b| ) or [Schneider et al.| ( |2011| l, 
who set the limit at around 25 Mjup, then 3 of the sub-stellar 
candidates are brown dwarf candidates (HD 1465 lb, HD30246b, 
and HD92320b), while the remaining two are planetary candi- 
dates (HD22781b and HD 156279b). A more detailed analysis of 
other possible classifications of the available objects is outside 
the scope of this paper and will be presented in a future work. 

The objects reported here increase the number of similar ob- 
jects present in the literature by about 12%, contributing to es- 
tablish solid obsei-vational evidence to test and constrain forma- 
tion and evolution theories. At the time of writing the number 
of sub-stellar objects or candidates more massive than 10 Mjup 
orbiting solar-type stars in relatively short iP < 10^ days) or- 
bits amounts to little less than 60 objects. The small size of this 
sample currently hinders most statistical analyses. Nevertheless, 
these objects warrant a rigorous statistical study, which might 
provide new and improved constraints for models and theories. 
Also, any effort to increase the number of objects in the mass 
range studied here, as well as all attempts to measure the actual 
mass of candidates, should be encouraged. 
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Table 2. Radial velocity measurements 



HD 14651 


BID 


RV 




BVS'' 


Exp. time 


S/N/pix 


-2 400 000 


(kms-i) 


(kms"') 


(km s~') 


(s) 


(at 550 rmi) 


54043.5099' 


CO 'O /I T 

52.243 


0.005 


-0.052 


500 


70 


^ A f\ A ^ A ^ ^ A \ 

54047.4534' 


51.611 


0.004 


-0.033 


300 


87 


54048.4953^ 


51.424 


0.004 


-0.044 


400 


94 


54049.4427^ 


51.220 


0.004 


-0.039 


500 


80 


54051.4619' 


50.745 


0.005 


-0.034 


400 


64 


54080.3977^ 


54.147 


0.004 


-0.033 


400 


69 


54501.2668 


53.781 


0.005 


-0.038 


227 


45 


54507.2945 


53.448 


0.005 


-0.040 


227 


46 


54705.6496 


53.075 


0.005 


-0.036 


341 


49 


54717.6278 


54.188 


0.005 


-0.032 


391 


47 


54720.6097 


54.228 


0.005 


-0.052 


364 


50 


55240.3093 


51.304 


0.005 


-0.028 


370 


48 


55405.5712 


49.679 


0.006 


-0.046 


227 


36 


55409.5920 


49.053 


0.005 


-0.030 


209 


53 


55453.6099 


53.805 


0.004 


-0.044 


188 


53 



HD16702 


BID 


RV 




BVS" 


Exp. time 


S/N/pix 


-2 400 000 


(kms-i) 


(kms-i) 


(kms-i) 


(s) 


(at 550 rmi) 


54365.6123 


4.844 


0.005 


0.005 


200 


54 


54407.4115 


5.699 


0.004 


0.003 


300 


79 


54435.4567 


5.275 


0.004 


0.006 


336 


81 


54461.4090 


2.086 


0.005 


0.015 


360 


44 


54462.3212 


2.091 


0.004 


0.009 


221 


60 


54463.3549 


2.141 


0.005 


0.018 


202 


53 


55066.6291 


6.403 


0.005 


-0.004 


279 


48 


55067.6359 


6.521 


0.005 


-0.001 


206 


51 


55068.6510 


6.663 


0.005 


0.003 


424 


47 


55069.6318 


6.742 


0.006 


-0.001 


600 


37 


55405.5953 


2.115 


0.005 


0.024 


347 


51 


55425.6437 


5.352 


0.005 


-0.008 


194 


47 


55427.6357 


5.811 


0.005 


0.004 


159 


50 


55436.6205 


6.896 


0.005 


-0.021 


285 


46 


55438.5971 


6.909 


0.005 


-0.010 


257 


46 


55444.6215 


6.572 


0.004 


0.004 


180 


62 


55449.5477 


6.089 


0.005 


-0.004 


202 


46 


55465.6213 


3.654 


0.006 


0.027 


270 


43 


55475.6539 


2.312 


0.005 


0.022 


189 


55 


55479.4840 


2.092 


0.005 


0.009 


276 


49 


55488.4233 


2.816 


0.006 


0.031 


200 


42 


55503.4724 


6.354 


0.005 


-0.009 


222 


47 
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Table 2. continued. 



HD22781 


BJD 


RV 


^ KV 


BVS° 


Exp. time 


S/N/pix 


-2 400 000 


(kms"') 


(kms"') 


(km s"^) 




(at 550 nm^ 

veil tJtJyJ 1.11.111 


54079.4169' 


8.154 


0.005 


-0.005 


600 


61 


54080.4125^ 


8.162 


0.004 


0.008 


500 


72 


54081.3843^ 


8.156 


0.004 


0.009 


500 


83 


54087.4543 


8.142 


0.005 


-0.005 


240 


43 


54097.4342 


8.134 


0.004 


0.011 


300 


56 


54125.3592 


8.096 


0.005 


0.001 


1000 


42 


54127.3144 


8.107 


0.005 


0.014 


680 


44 


54133.2899 


8.085 


0.004 


-0.004 


241 


60 


54138.3552 


8.095 


0.004 


0.003 


724 


60 


54141.3158 


8.085 


0.004 


0.009 


483 


75 


54148.3446 


8.083 


0.005 


0.021 


241 


41 


54152.3596 


8.086 


0.006 


-0.003 


241 


48 


54173.3081 


8.063 


0.004 


0.007 


463 


67 


54344.6620 


8.297 


0.004 


0.014 


367 


64 


54372.6427 


8.994 


0.005 


0.010 


700 


40 


54375.6321 


8.905 


0.003 


0.009 


734 


97 


54376.6144 


8.899 


0.003 


0.002 


734 


102 


54381.4955 


8.816 


0.005 


0.009 


194 


47 


54406.4668 


8.604 


0.004 


-0.002 


600 


68 


54407.5038 


8.571 


0.004 


0.002 


400 


74 


54409.5721 


8.577 


0.004 


0.007 


360 


59 


54855.3067 


7.990 


0.005 


0.009 


831 


47 






n oos 


014 




4Q 


54877.2718 


8.715 


0.005 


0.008 


558 


48 


54878.3020 


8.831 


0.005 


0.005 


638 


49 


54890.2711 


9.297 


0.004 


0.017 


288 


52 


54893.3243 


9.211 


0.005 


0.001 


696 


46 


55079.6628 


8.248 


0.007 


0.015 


900 


33 


55088.6638 


8.239 


0.005 


0.009 


552 


47 


55214.3938 


8.101 


0.005 


0.006 


350 


49 


55409.5967 


9.222 


0.004 


0.006 


291 


52 


55444.6436 


8.752 


0.004 


-0.010 


218 


53 
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Table 2. continued. 



HD30246 


BJD 


RV 


O-RV 


BVS" 


Exp. time 


S/N/pix 


-2 400 000 


(kms-i) 


(kms-i) 


(km s-i) 


(s) 


(at 550 nm) 


54051.5549' 


41.300 


0.006 


0.020 


400 


61 


54053.5 157^ 


41.301 


0.004 


-0.002 


350 


78 


54078.4650^ 


41.231 


0.007 


-0.011 


500 


53 


54079.4428'" 


41.240 


0.005 


-0.009 


500 


67 


5408 1 4249'" 


41.250 


0.004 


-0.009 


400 


91 


54088 4622 1^ 


41.234 


009 


-0 029 


233 


44 




41 1 Q8 


U.UUJ 


-U.UIZ 




D / 




/111 ^'2 


U.UU/ 


n (\r\A 
U.UU4 




J 1 


54127.3579' 


41.167 


0.006 


-0.022 


680 


57 


54133.3531' 


41.108 


0.006 


0.032 


180 


58 


54148. 3614"^ 


41.067 


0.006 


0.004 


180 


53 


54152.3558^ 


41.024 


0.008 


-0.005 


180 


54 


54173.3335^ 


41.005 


0.005 


0.010 


295 


74 


54855.3659 


41.668 


0.009 


-0.043 


1264 


30 


55096.6444 


41.180 


0.008 


0.005 


295 


33 


55214.4034 


40.810 


0.005 


-0.015 


216 


48 


55240.3314 


40.733 


0.006 


0.015 


355 


47 


55268.2826 


40.585 


0.005 


0.010 


209 


50 


55283.3259 


40.592 


0.007 


0.024 


236 


39 


55289.3131 


40.635 


0.005 


0.002 


367 


60 


55437.6466 


42.800 


0.006 


0.007 


603 


46 


55465.6651 


42.684 


0.005 


-0.017 


376 


51 


55488.4488 


42.570 


0.006 


0.003 


474 


50 



HD92320 



BJD 


RV 


O-RV 


BVS" 


Exp. time 


S/N/pix 


-2 400 000 


(kms-i) 


(kms-1) 


(km s-i) 


(s) 


(at 550 nm) 


54090.6482 


1.911 


0.004 


-0.022 


300 


66 


54093.7075 


1.542 


0.004 


-0.016 


300 


71 


54096.6818 


1.174 


0.004 


-0.015 


300 


64 


54097.6492 


1.071 


0.005 


-0.013 


300 


52 


54584.4482 


-1.845 


0.005 


-0.014 


180 


44 


54601.3449 


-1.957 


0.005 


-0.024 


885 


51 


54642.3592 


1.363 


0.005 


-0.011 


309 


50 


54645.3648 


1.939 


0.005 


-0.028 


295 


47 


54889.4636 


-1.933 


0.005 


-0.007 


638 


46 


54926.4396 


0.218 


0.005 


-0.020 


192 


50 


54927.5196 


0.387 


0.005 


-0.011 


236 


49 


55266.5205 


0.467 


0.005 


0.003 


288 


47 


55305.3537 


-1.718 


0.005 


-0.017 


257 


51 


55519.7222 


2.247 


0.005 


-0.009 


185 


48 


55525.6987 


3.046 


0.005 


-0.016 


365 


47 


55527.7043 


3.184 


0.005 


-0.016 


270 


45 
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Table 2. continued. 



HD137510 




KV 


O'RV 


£> V O 


Exp. time 


C /NT /triv 

O/iN/plX 


O A C\f\ f\C\C\ 

-I 4UU UUU 


(Km S ) 


(Km s ) 


(km s ) 


(s) 


(at jjU nm) 


54549.6357 


-6.711 


0.006 


-0.019 


245 


57 


54553.5731 


-6.707 


0.006 


-0.003 


128 


56 


54554.5767 


-6.709 


0.004 


-0.015 


180 


100 


54873.6698 


-6.023 


0.005 


0.008 


192 


79 


55051.3403 


-6.285 


0.007 


-0.046 


41 


50 


55051.3426 


-6.267 


0.005 


-0.009 


181 


80 


55240.7203 


-6.753 


0.007 


-0.035 


35 


48 


55266.7010 


-6.754 


0.007 


-0.003 


65 


47 


55268.6309 


-6.777 


0.007 


0.012 


37 


49 


55288.5197 


-6.766 


0.004 


-0.013 


194 


102 


55305.3824 


-6.756 


0.005 


-0.004 


180 


83 


55391.3744 


-6.699 


0.006 


-0.012 


35 


53 


55436.3653 


-6.637 


0.006 


-0.025 


180 


61 



HD 156279 


BJD 


RV 


o-Rv 


BVS" 


Exp. time 


S/N/pix 


-2 400 000 


(kms-i) 


(kms-i) 


(kms-i) 


(s) 


(at 550 nm) 


55330.5522 


-20.627 


0.004 


-0.036 


327 


48 


55361.4611 


-20.881 


0.004 


-0.031 


148 


50 


55398.3821 


-20.163 


0.003 


-0.037 


377 


90 


55405.3572 


-20.142 


0.004 


-0.032 


191 


50 


55424.3514 


-20.370 


0.004 


-0.035 


114 


51 


55425.3640 


-20.360 


0.004 


-0.017 


250 


48 


55427.3726 


-20.376 


0.004 


-0.028 


176 


48 


55441.4021 


-20.502 


0.004 


-0.035 


173 


51 


55449.2932 


-20.541 


0.004 


-0.029 


193 


48 


55479.2757 


-20.766 


0.005 


-0.042 


180 


45 


55503.2440 


-20.999 


0.004 


-0.042 


182 


49 


55513.2168 


-21.137 


0.004 


-0.035 


221 


65 


55525.2393 


-20.823 


0.005 


-0.042 


502 


45 


55540.2372 


-20.210 


0.005 


-0.025 


806 


48 


55584.7229 


-20.592 


0.004 


-0.031 


180 


54 



HD167215 



BJD 


RV 




BVS" 


Exp. time 


S/N/pix 


-2 400 000 


(kms-i) 


(kms-1) 


(km s-i) 


(s) 


(at 550 nm) 


55042.3613 


-44.609 


0.008 


0.075 


207 


49 


55051.3730 


-44.845 


0.008 


0.068 


198 


48 


55054.3550 


-44.932 


0.008 


0.028 


181 


50 


55067.4087 


-45.354 


0.008 


0.021 


195 


48 


55331.6227 


-45.066 


0.008 


0.047 


209 


46 


55333.5844 


-45.022 


0.008 


0.055 


408 


45 


55366.5956 


-44.676 


0.010 


0.051 


192 


39 


55378.5849 


-44.582 


0.008 


0.039 


198 


49 


55409.4168 


-44.340 


0.008 


0.051 


183 


49 


55423.3895 


-44.211 


0.008 


0.040 


425 


45 


55425.3973 


-44.200 


0.008 


0.065 


161 


46 


55429.3835 


-44.176 


0.008 


0.022 


114 


49 


55448.3524 


-44.062 


0.008 


0.003 


209 


47 


55479.3607 


-43.890 


0.010 


0.013 


225 


44 



Bisector Velocity Span 

Observations in thosimult mode before upgrade (see 5 2.2 1. 



